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An exact solution for an XXZ chain with four-body interactions is obtained and its phase dia- 
gram is determined. The model can be reduced to two chains coupled by four-body interactions, 
and it is shown that the ground state of the two-chain model is magnetized in part. Furthermore, 
a twisted four-body correlation function of the anti-ferromagnetic Heisenberg chain is obtained. 
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§1. Introduction 

Quantum spin systems solved exactly so far are es- 
sentially one dimensional. It is a challenging problem 
to find a two, or higher, dimensional exactly solvable 
quantum spin model. As a first step we study exactly 
solvable models for two, or more, chains combined each 
other by some interaction^. A well-known example is the 
Majumdar- Ghosh model Id-* which can be interpreted as 
a two-leg zigzag ladder; 



H, 



MG 



1 ^ 1 ^ 



N 



J = l 



N 



+ T.^rT,+J2TrS,+u (1.1) 

where Sj and Tj are operators of spin- 2 at site j on each 
chain. The ground states are given by doubly degenerate 
dimer states with no magnetization. 

We introduce here a novel two-chain model, which is 
exactly solvable by the Bethe ansatz method. The model 
consists of two Heisenberg chains coupled by four-body 
interactions. It can be expressed as twisted two chains; 

1 ^ 1 " 

i?2chain ^ jZ^^J ' -^i+l + 7 Z^ ^J ' ^J + 1 



N 



J = l 



^(5,xT,)-(5,+ixT,+i) 



(1.2) 



N 
J = l 

This Hamiltonian is obtained as a special case of the 
following single-chain model; 

N N 

N 

+ Js^yiSj ■ Sj+2){Sj+i ■ Sj+3) 



- [Sj ■ Sj+3]a3[Sj+i ■ Sj+2]a'^ 

+ const., (1-3) 

where an abbreviation of anisotropic scalar product is 
used, 



[a ■ b]A = a^6^ + aHy + Aa'¥ 



(1.4) 



This model obtains integrability when the anisotropics, 
Ai, A3, and A3, and the strengths of interactions, Jk, 
satisfy certain conditions with two parameters, A and a. 
In the next section the integrability condition of the 
single-chain model (1.3) is presented. The model is in- 
tegrable if its Hamiltonian is decomposed to a sum of 
mutually commuting operators including the XXZ model 
Hamiltonian. The eigenvalues of the commuting opera- 
tors are evaluated in §3. The existence of phase tran- 
sition is discussed in §4, where the phase boundary is 
determined analytically. The magnetic property of the 
two-chain model is mentioned briefly in §5. In §6 a corre- 
lation function of the nearest four sites of the Heisenberg 
chain is explicitly given. The last section is devoted to a 
brief summary. 

§2. Integrability 

We consider the system of the following Hamiltonian; 

HA,a ^ Hxxz + aHint, (2.1) 

where the first term on the r. h. s. is the Hamiltonian of 
the usual anti-ferromagnetic XXZ model, 

^ A 

Hxxz - J2^S, ■ S,+i]a - -jN, (2.2) 

and the interaction term, iJint, involves up to the third 
neighbor interactions; 



1 + A^ 



J = l 



2A 

Ta^' 



N 



~^Z^^J' ^i- 



j=i 
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N 

- [Sj ■ S,+3]a[S,+i ■ 5,+2]a-i) + jN. (2.3) 

The coupling constant, a, takes a real value, the 
anisotropy parameter is restricted in the range A > — 1, 
and the periodic boundary condition, Sjy-^-i — Si, i s im- 
posed. It is easy to see that the Hamiltonian (2J.) has 
a similar form to that of the single-chain model (1.3). 
The former can be considered to be a special case of the 
latter. 

Note that the parameter a does not denote the 
strength ratio of the second-neighbor, or four-body, in- 
teraction to the nearest neighbor one, but expresses 
only the ratio between Hi^t and Hxxz- The strengths, 
Jfe, and the anisotropics, Aj,, of the interactions in the 
Hamiltonian (1.3) should be related with the two param- 
eters, a and A, as follow; 



Ji = l 



1 + A2 



Aa 
~2 



a, J2 = -^, J3 = 2Ac 



Ai 



2A(l-a) 
2-(l + A2), 



A., 



(a:,)- 



(2.4) 
(2.5) 



The strength ratio of the second-neighbor interaction to 
the four-body one, J2/J3, is fixed at 1/4, though J2/J1 
has freedom. 

Conversely, the Hamiltonian (1.3) is integrable un- 
der these conditions, (2.4) and (2.5), since it is decom- 
posed to a sum of two commuting Hamiltonians, ( |2.2| ) 
and (2.3). The reason of the commutativity is that the 
complicated interaction operator ifi„t is introducecLby 
the help of the quantum inverse scattering methocH'H^ 
as follows. 

The Hamiltonian of the XXZ model is related with 
the transfer matrix, r(A), of the six- vertex model. The 
transfer matrix is given as the trace of the monodromy 
matrix; 

r(A)=trr(A). (2.6) 

The monodromy matrix is constructed of L-matrices; 

r(A) = LAr(A)iN_i(A)---Li(A). (2.7) 

The L-matrices are defined as follow; 

1 +0'aiix<8crf 1 - CTa„„(g)(7f 

L,W - «(A) '^^ + d{X) ^^^ 



sinh^(A + i) 



tt(-: 



aux-JCTj +craux«)Cr|), (2.; 



with 



a(A) 



1 



d{X) = -7 



1 sm 7 

1 sinh ^(A — i) 



isin7 sinh ^(A -I- i) ' 



(2.9) 



(2.10) 



where the parameter 7 takes a real or a purely imaginary 
value, and a"^ = \{(t^ ± cr^). The subscript of the Pauli 
matrices, cr — {a^, a"^ ,a^), denotes the space upon which 
they act. The trace and the products of L-matrices in 
the definition of the transfer matrix should be performed 



in the auxiliary space of the operator cTaux- 

The L-matrices satisfy the Yang-Baxter relation; 

i?(A,/i)(L(A)(»L(Ai)) = (i(/i)®L(A))i?(A,M), (2.11) 
where the _R-matrix is given as 



i?(A,/i)=/(A,A*)- 



gi^fj-)- 



with 



/(A,/i) = 
5(A,/i) 



sinh ^ (A — /i + 2 i) 
sinh ^(A — /i) 
isin7 



(2.12) 

(2.13) 

(2.14) 



sinh ^ (A — /i) 

Then the transfer matrix forms a commuting family with 
the spectral parameter A; 



[r(A),r(/i)]=0. 



(2.15) 



The fc-th order Hamiltonian, _ff ^'^^ , is defined as a log- 
arithmic derivative of the transfer matrix; 



hW 



sin 7 d 
~^dX 



InT(A) 



(2.16) 



The XXZ model Hamiltonian ( p.2[ ) is identical to the 
Hamiltonian of the first order; 



H 



XXZ 



= H^^\ A = cos7. 



(2.17) 



Because of the commutativity of the transfer matrix, 
the Hamiltonians of arbitrary order also commute with 
one another; [H'^^\ H'^^'^] = 0. Then they have common 
eigenstates, which are especially common to those of the 
XXZ model. A linear combination of the commuting 
Hamiltonians, 



H^J2^>^H^'^^ 



(2.18) 



fc=i 



is, therefcwje, integrable by the Bethe ansatz method. 

TsvelikcP and FrahnJa^ investigated a model con- 
structed of the first and the second order Hamiltonian; 



H :^H 



XXZ 



aH^^\ 



(2.19) 



where the explicit form of the interaction operator is 
given as 

N 

i?(2) = A ^[5, • (5j_ix5j+i)]a-i. (2.20) 

They showed that this model exhibits a phase transition 
at some finite value of the coupling constant. The crit- 
ical value, acrit, for arbitrary anisotropy is analytically 
obtained by Frahm. The ground state at \a\ > acrit ob- 
tains an incommensurability and it is also magnetized in 
the Ising-like region; A > 1. 

The model, however, has the strange interaction (2.20) 
which breaks the parity invariance. That motivated us 
to choose the new interaction (2.3), which is nothing but 
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the third order Hamiltonian; 



for A: > 1; 



(2.21) 



This operator conserves the parity. Furthermore, it con- 
tains isotropic second- neighbor two-body interactions. 
When we take the parameter as a = 1 and A = 1 in 



Ekii^,}) 



sin 7 d 
7 dn 



fe-i 



E{^,,{x,}) 



with 



the Hamiltonian (2.1), the nearest neighbor interactions 
cancel out. It becomes the following form; 



S(/i,{A,}) = ^|-ln0(M,{A,}). 



(3.5) 



(3.6) 



H,i^ 



1 ^ 



'j+2 



N 
l^(5,xS, + i).(5,+2X5,+3) 
J = l 



N 



(2.22) 



The second term on the r. h. s. of the eq. ( |3.4D docs not 
contribute to the eigenvalue of the Hamiltonian if fc < iV, 
since the function {d{^))^ and its derivatives vanish at 
fj, = i. We consider the thermodynamic limit, iV — > oo, 
and the term can be omitted. Then the function ( |3.6| ) 
reduces to 



This is equivalent with the Hamiltonian of the two-chain 
model; Hi,i = 2iJ2chains - N/8. 

§3. Eigenstates and Eigenvalues 

In this section we investigate the system of the 
by using the algebraic Bethe 



M 






with 




the 
and the thermodynamic Bethe ansatz 



Hamiltonian 
ansatz rp.e± 
method.ErlZl' 

We are interested in the ground state of the system. 
Then we consider the low temperature limit to observe 
the ground state. 

Since the system has the same eigenstates as those of 
the XXZ model, the states with S^^^^^ = (iV-2Af)/2 are 
determined by the Bethe ansatz equations of the XXZ 
model. Using the algebraic Bethe ansatz method, we 
obtain the following equations; 



e(^,A) 



sin^ 7 



(3.7) 



(3.8) 



cos 7 — cosh7(A ~ n + i) 

The eigenvalue of the fc-th order Hamiltonian is given as 
the derivative of the function (3.7); 

M 

^fc({A,})=E'^^-(^^)' (3.9) 



i=i 



with 



efc(A) 



sin 7 d 



fc-i 



e(i,A). 



(3.10) 



a(A,) 
d{X,) 



N 



n 



/(A,-,AO 
/(A/,A,r 



7 dX, 
Note that the following identity is used in the derivation 



j = l,2,...,M. (3.1) of thecq. (|3.iq) 



where the functions a(A), d{X), and /(A,/i) are given in 
the eqs. (2.9), (2.10) and (2.13). The functions on the 
1. h. s., (a(A))''^ and (d(A))^, are the eigenvalues of the 
diagonal elements of the monodromy matrix, [T(A)]i.i 
and [r(A)]2,2, respectively, corresponding to the pseu- 
dovacuum state; | ) = Ylj^i I T )i- 

An eigenstate corresponds to a set of the spectral pa- 



^^ei^i,Xj) = ~-—e{^,Xj] 



d^j. 



dX, 



(3.11) 



The expression (3.£) means that the system is com- 
posed of M particles, each of which has rapidity A, and 
the energy of the system is given as the sum of one- 
particle energies (3.10). 



§4. Phase Diagram of the Single-Chain Model 



rameters, {Aj}, which is a solution of the eqs. (3.1). The 



state vector can be expressed by using the operator of 
an off-diagonal element of the monodromy matrix; 



We determine the phase diagram of the system (2.1) 
here in the thermodynamic limit. The energy for the 



M 



i{A,})=n[^(^^)]i.2i 



eigenstates (3.2) is given from the results of the previous 
section; 



(3.2) 



M 



Now we calculate the eigenvalue of the fc-th order 
Hamiltonian ( p. 16 ) corresponding to the state (3.2). The 
eigenvalue of the transfer matrix is calculated as follows; 



^ = E 



/ sinh 7 d 
[~^dX, 



sin 7 



cos 7 — cosh 7Aj 



(4.1) 



T(M)|{Aj})=e(M,{A,})|{A,}), 



(3.3) 



where 

M M 

0(/i, {A,}) = (a(^))~ n /(^^■' m) + km))^ n /(a^' ^:'-)- 

(3.4) 
Taking a logarithmic derivative of this equation, wc ob- 
tain the eigenvalue of the k-th order Hamiltonian H'-''^ 



The parameter 7, which is related with the anisotropy as 
A = cos 7, takes a real value when |A| < 1 and a purely 
imaginary value when A > 1 . 

4.1 |A|<1 

Here we consider the case that the absolute value of 
the anisotropy is less than one. The Bethe ansatz equa- 
tions (3.1) have string solutionaa'B'' in the thermody- 
namic limit N —^ 00; 



X) = Xj + {uj -2k + l)i + 



1 



27 



-TTl, 



Norihiro MuRAMOTO and Minoru Takahashi 



for fc = 1, . 



(4.2) 



where the length, rij, of a string is bounded by the num- 

of a 



ber of the particles, X^^-j = ^^■ 

string takes value of ±1. The parities originate from the 



and the parity, Vj, 



periodicity of the eqs. (3.1) along the imaginary axis. It 
is thought that the string solutions are sufficient to form 
thermodynamics. _ j.. 

Following the usual procedure, B'cP we classify the so- 
lutions of the Bethe ansatz equations by the lengths and 
the parities of strings. The density of rapidity and that 
of holes, Pj{X) and PjW, where the subscript denotes 
a class, is defined. In the thermodynamic Bethe ansatz 
method, the dressed energy, Sj{X), is introduced as 



e,(A)=Tlnr/,(A), rj,{X) 



pMY 



(4.3) 



We can deduce a set of functional equations about pj (A) , 
p^(A), and 77^ (A) under the thermodynamic equilibrium 
condition, but it is lengthy and cumbersome to write 
down them and they are omitted here. (See refs. |g and 
g if necessary.) 

Though the number of the classes is infinite, in the low 
temperature limit we only have to consider two of them, 
both of which belong to the length one, one belongs to 
the positive parity and the other to_thc negative, since 
it is shown by Takahashi and SuzukiB' that the energies 
of the classes belonging to the lengths longer than one 
are large and that they don't contribute to the ground 
state. Then we obtain the following functional equations 
about the densities of rapidity, p±(A) and pij.(A), and the 
dressed energies, e±(A), of the two classes; 

g±{\) = ±{p±{X) + pl{\)) 

+ K±* p+iX) + K^* p^{\), (4.4) 

e±(A) = 4(A) - K±* 4-) (A) + K^* e^_r\x), (4.5) 



with 



5±(A) = ±- 
K±{X)=± 

4(A) 



7Sin7 



27r cosh 7 A =F cos 7 
1 7 sin 27 

27r cosh 7A ^ cos 27 ' 



2tt sin 7 






(A) 




1 



sin 7 d 



7 

£±(A) < 0, 
otherwise. 



dX 



5±(A), 



and the asterisk product a * b{X) denotes a convolution 
defined as 



a * b{X) 



a{X — p) h(p) dp. 



(4.6) 



Note that in the low temperature limit, T ^ 0, the fol- 
lowing relation holds from the definition ( }4.3| ): the ra- 
pidity density p{X) (the hole density p^{X)) is zero when 
the dressed energy takes a positive (negative) value. The 
eqs. (4.4) together with this relation become complete to 
determine the densities of rapidity. 

The analytical solutions of the dressed energies are 



easily obtained if the coupling constant \a\ is sufficiently 
small. By applying the Fourier transformation to the 
eqs. (4.5) we obtain 



e±(A) 



7rsin7 



1 



sin 7 d 
7 dX 



/±(A), (4.7) 



with 
/+(A) =sech— , 







<7< 



/-(A) = 



cos(7r-^/27) cosh(7rA/2) n 



- < 7 < TT. 



coshTrA + cos{'K^/j) 

The solutions ( [l.TJ ) are valid as long as e+(A) < and 
e_(A) > for arbitrary value of the rapidity A. (Note 
that — cos(7r^/27) is positive when tt/2 < 7 < tt.) It is 
easy to see from the solutions that this condition is sat- 
isfied when the coupling constant lies in the finite range, 
"^rit < " < "iif where 



a: 



with 



h+{x) = 



h (x) 



crit 

1 
1 



± 



27 



TT sin7 



KiTr'h), 



(4.8) 



0<7<2- 



sin^(a::/2) tt 



l-|-sin^(a;/2) 



2 sin^ X 



2 < 7 < Ti", 



Btt 



<7 < 

-— < 7 < TT. 

5 



2 cos^ X — 4 cos X + 3 

The ground state is the same as that of the XXZ model 
when the coupling constant is in this range, since the 
Fermi points of the dressed energies are the same. When 
the coupling constant goes beyond the range, there add 
new four Fermi points. 

The appearance of the new Fermi points is easily ob- 
served at the 'free fermion point', A = 0, where the XXZ 
model reduces to a free fermion system. The one-particle 
dispersion of our model is easily found to be 

e(p) = (1- ■^)cosp- -cos3p. (4.9) 

There appear three negative-energy regions with six 
Fermi points when the coupling constant is beyond the 
critical value, \a\ > a^.j^ — 1. 

Similar situation can be observed on the two dressed 
energies at another value of the anisotropy parameter. 
Then we conclude that there occurs a phase transition if 
the coupling constant goes across the critical value (p~^). 
The new phase is characterized by added new branches 
of excitations, which are also gapless. 

The ground state in the new phase has no magnetic 
moment. We can see this fact by taking a Fourier com- 
ponent of one of the eqs. (4.4) at zero frequency. It leads 
to 

M/N = S^,,JN 
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= 1- J P+i\)d\- I p.iX)dX 



,^jpHx)dx o<,<-, 



(4.10) 



Since it can be proved that e_i^\{X) < everywhere for 
< 7 < 7r/2 {tt/2 < 7 < tt), there is no hole of rapidity 
belonging to — (+)-parity. The ground state is, therefore, 
not magnetized at any strength of the interaction (2.3). 
The new phase is similar to the corresponding one of the 
Tsvelik-Frahm modeHQ' in this point. 

4.2 A>1 

Here the parameter 7 is purely imaginary and we put 
J = i9, 6 > 0. Since in this case there is no periodic- 
ity along the imaginary axis in the Bethe ansatz equa- 



tions (3.1), only string solutions with positive parity, 
which are centered on the real axis, are allowed. On 
the other hand, the equations become periodic along the 
real axis. We restrict, therefore, the rapidities in the 
following range; 



A: 



{A|-^<A<^}. 



(4.11) 



(When the parameter 7 is zero, an adequate limit, — > 0, 
should be taken on the equations below in this section 
and the range of the rapidity should be expanded to in- 
finity; A = {A|-oo<A< 00}.) 

We only have to consider string solutions belonging to 
the length one in the same way as the previous section. 
The density of rapidity, p(A), and the d: 
e(A), are determined by the following eqs 



g{X)^p{X)+p\X)+K*p{X), 
e(A)=e°(A)-/^*e(-)(A), 



energy, 

(4.12) 
(4.13) 



with 



5(A)- 
K{X) = 



61 sinh 6* 



2tt cosh 9 — cos 9X ' 
1 9 sinh 29 



s"(A) = - 



^(-)(\) = 



27rcosh26' 
27r sinh 9 



- cos9X 

/sinh0 d 
1 



(A) 




£(A) < 0, 
otherwise, 



dX 



2n 



5(A), 



and the asterisk product a * 6(A) denotes a convolution 
over the range (4.11) here; 



a ^ b{X) — / a{X — p)b{p) dp. 
Ja 



(4.14) 



The eq. (4.13) has an analytical solution for sufficiently 
small |a|. The dressed energy is obtained as 



e(A) = - 



1 + a 



/sinh 61 d 



dX 



1 



Yl l^sech\n\9 



K{k)smh9 



xdn{K'{k)X;k), 



/sinh 61 d 
\ 9 dX 



(4.15) 



where the parameters K{k) and K'(k) are the complete 
elliptic integrals of the first kind with modulus k and 
k' = \/l — fc^, respectively, and the function dn(x; k) is 
one of Jacobi's elliptic function with modulus k. The 
modulus is related with the anisotropy as 

_K'{k) 



cosh0, 9 = 



K{k) 



(4.16) 



The solution (4.15) is valid as long as £(A) < ev- 
erywhere. This condition is identical with the following 
one; 

The critical value is continuous at 7 = with that of the 
previous section taking value of ±(2/7r)^. 

The ground state coinci des w ith that of the XXZ 
model when the condition (4.17) is satisfied. It is not 
magnetized. If the anisotropy is in the Tsing-like region', 
A > 1, the ground state also shows the Neel order and 
has excitations with gap. 

When the coupling constant is beyond the critical 
value, the dressed energy gains a positive region where 
e(A) > 0. We name the region 13+ , and define its com- 
plimentary set D^; 



D+ = {Ae A|e(A) >0}, D'^K-D^ 



(4.18) 



The existence of not-empty region D^ means that there 
appear new Fermi points and the system has gapless ex- 
citations, since the dressed energy function, e(A), is con- 
tinuous and takes a definite value at the end points of 
domain A. Thus the system shows a phase transition 
at the critical value of the coupling constant. When the 
anisotropy parameter lies in the Tsing-like region', the 
system undergoes a transition from the gapful phase to 
the gapless one when the coupling constant exceeds the 
critical value. 

Another typical character of the new phase is that the 
ground state is magnetized. We can calculate the mag- 
netic moment as 



MIN = Sl,JN 

4-/^p(A)dA 

ph(A)dA, 



(4.19) 



by taking the zero frequency component of the Fourier 
expansion of the eq. (4.12). As we noted in the previous 
section, the hole density, /o'^(A), can have non zero value 
only in the region £)+. When the coupling constant is 
large enough for the region D"*" not to be empty, it is easy 
to make sure from the eq. (4.12) that there exist holes. 
Thus the last term of the eq. (4.19) gives a positive value 
and the system obtains a magnetic moment. 

Note that the system is not completely but partly mag- 
netized even for quite large coupling constant, since the 
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region D never becomes empty. In fact the magnetic 
moment gets saturated at S^^^^JN = 0.3217 (0.1479) 
when a -^ +00 (—00) and A = 1. 




anisotropy: A 



Fig. 1. Phase diagram of the single-chain model (2.1) is given. Its 
ground state has the same property as that of the XXZ model 
for sufficiently small |o| between «^. and Oi~ . ; the phase (I), 
A < 1, has gapless excitations, the phase (II), A > 1, has exci- 
tations with gap and the ground states in both phases are not 
magnetized. The systems in the region (III) and (IV) obtain 
additional branches of excitation to those in phase (I) and (II), 
respectively. The extra excitations in both phases are gapless 
ones. The ground state in the phase (III) is also not magnetized, 
though in phase (IV), A > 1, it obtains a finite, but not-full, 
magnetization. When the parameters are chosen as A = q = 1, 
the model reduces to the two-chain model, which lies in the par- 
tially magnetized phase (IV). 



§5. Two-Chain Model 

Now we can get information of the two-chain model. 
We only have to watch the single-chain model (2.1) at a 
special value of the parameter, a ~ 1 and A = 1. From 
the view point of the single-chain model, the coupling 
constant is beyond the critical value, acrit = (2/7r)^ < 1. 
Then the system of the two-chain model stays in the 
partly magnetized phase. 

The magnetic moment is calculated numerically; 



^total 



/N = 0.1049, 



(5.1) 



which shows that the system is about 20% magnetized 
in the ground state. As noted in the previous section, 
the system has gapless excitations with extra branches. 

§6. Twisted Four-Body Correlation Fiinction 

The expectation value of the fc-th order Hamiltonian 
operator, if ^'^^ for the ground state of the isotropic and 
anti-ferromagnetic Heisenberg model {XXX model) is 
presented in this section. We evaluate it by using the 
eq. ( ^.9| ) and the rapidity distribution function, po(A), of 
the ground state of the XXX model; 

po(A) = ^sech— 



1 f°° 
— / sechcL;exp(— iAtj) do;. 

47r J_oo 



(6.1) 



Note that in this section we analyze only the isotropic 
case, A = 1. 

The expectation value, (if ''^'-')o, is calculated from the 
eq. ( |3.9[ ) in an analytical form; 

(ii«)o 



E, 



N 



(6.2) 



This vanishes if the integer k is even, and it gives the 
following value for odd k; 

/oo 
efe(A)po(A)dA 
-00 

(_l)(fe+i)/22-fe+i(i_2-'^'+i) 
x(fc-l)!C(fc) forfc>3, (6.3) 

-In 2 forfc = l, 

where the function ^(s) is Riemann's zeta, C(s) = 

The expectation value of the first order Hamiltonian 
leads the nearest neighbor correlation function of the 
XXX model; 



{Sj ■ Sj+i)a - 



(if(") 

N 



1 



= i-ln2. 
4 



(6.4) 



Similarly the expectation value of the two-chain model 
Hamiltonian (2.22) gives the following relation; 

-{Sj ■ Sj+2)o + '2{{Sj X Sj+i) ■ {Sj+2 X Sj+3))o 

= i-ln2 + ^C(3). (6.5) 

The correlation function between the^econd-neighbor 
sites is obtained by one of the authors ;E1P 

(5,.S,+2)o = J-41n2+^C(3). (6.6) 

Then the twisted four-body correlation function is cal- 
culated as 

{{S,xS,+,) ■ (5,+2x5,+3))o = iln2- ^C(3) 

= -0.104198, (6.7) 

or expressed by its component as 

Wj •^j+l^j+2^j+3 ^j ^] + l'^]+2'='j+3/0 

^C(3)-^hi2^(l-<5„^) 

= 0.0173663 (l-(5c«/3), a,(3^x,y,z. (6.8) 

To check the validity of this result, we calculate the four- 
body correlation function for finite systems using the ex- 
act diagonalization method and plot them. Figure. 2 
supports that the eq. (6.8) is correct. 



§7. Summary 

The XXZ model with competing four-body interac- 
tion has been studied. When the four-body interaction 
is weaker than the critical value, acrit: the ground state 
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Fig. 2. The twisted four-body correlation functions of the anti- 
ferromagnetic Heisenberg chain with N sites, A'^ = 10, 12, 14, and 
16, are plotted as a function of l/N"^. In the thermodynamic 
limit, N — » oo, these values approach the theoretical value (6.8). 



coincides with that of the XXZ model and the spectrum 
of low-lying excitation is also the same. If it gets stronger 
than Qfcrit, the new phases arise. The critical value of the 
coupling constant has been obtained analytically. The 
new phases (III) and (IV), see Fig. 1, have gapless exci- 
tations. The phase (III), |A| < 1, the ground state is not 
magnetized and the phase (IV), A > 1, it is magnetized 
in part. 

The two-chain model (1.2) coupled by four-body in- 
teractions can be deduced from the above single-chain 
model. Since the former model corresponds to the lat- 
ter in the phase (IV) , the ground state of the two-chain 
model shows the partly magnetized property and has ex- 
citations without gap. 

From the expectation value of the four-body interac- 
tion operator, the twisted four-body correlation of the 
anti-ferromagnetic Heisenberg model has also been ob- 
tained analytically. 
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